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THEORETICAL INVESTIGATION OF A PROPORTIONAL-PLUS -FLICKER 

AUTOMATIC PILOT^ 

By Ernest C. Seaberg 


SUMMARY 


The proportional -plus-flicker automatic pilot operates by a non- 
linear principle whereby a fast -acting flicker servomotor response is 
combined with a low-speed proportional servomotor response for the pur- 
pose of obtaining supersonic stability and control. Essentially ; the 
autopilot maintains a zero reference about which the output is propor- 
tional to the input. However^ a flicker response overrides this propor- 
tional response at a fixed angle of gimbal displacement on either side of 
the zero gyroscope reference. Therefore^ in contrast to other high- 
speed control systems, the design requirements are simplified because 
the two components of the proportional-flicker control system are easy 
to build separately and they can be combined in a relatively simple 
manner. 

By application of the proportional-flicker principle, satisfactory 
stability can be obtained by the proper adjustment of the variable fac- 
tors in the autopilot mechanism, namely, the proportional gain, the 
amplitude of flicker -control deflection, the autopilot time-lag factor 
(the time lag between flicker and proportional operation), and the point 
in the range where the autopilot switches from a flicker to a proportional 
system. There is a possibility that these factors can be adjusted so 
that a more rapid response time (the time to reach steady state) is 
obtained with the nonlinear proportional -flicker autopilot than with a 
purely linear proportional autopilot. 

For the main part of this analysis, the proportional part of the 
system is approximated by a zero-phase-lag proportional autopilot with 
the assumption that the control surface moves instantaneously at the 
point where the system switches from flicker to proportional. Good cor- 
relation is shown between the results obtained by this method and results 
obtained by using a close approximation of an actual autopilot transfer 
function for proportional autopilot operation. 



^Supersedes recently declassified NACA RM L50I19^ 1950* 
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The proportional-flicker control system appears to he a practical 
method for obtaining pitch stabilization of a supersonic pilotless air- 
craft. Therefore^ trials of this system, particularly in supersonic 
vehicles, appear warranted. 


INTRODUCTION 


As part of the general research program for investigating various 
means of automatic stabilization, the Pilotless Aircraft Research 
Division of the Langley Aeronautical Laboratory has been conducting a 
theoretical analysis to determine the feasibility of using a proportional- 
plus -flicker automatic pilot for stabilization and control of a super- 
sonic canard airframe. The principle of operation of this autopilot is 
believed to be unique because it combines a fast -acting flicker servo 
response with a low- speed proportional servo response in a relatively 
simple manner, that is, by overriding the proportional part of the 
system at a fixed angle of gimbal displacement from the zero gyroscope 
reference through the use of simple electrical pickoffs attached to the 
displacement gyroscope outer gimbal. 

The reason for attempting to develop an autopilot of this type is 
to overcome the apparent difficulty in building a high-speed proportional 
servomotor by the construction of a proportional -flicker servo. Because 
the fast -acting flicker servo and the low-gain, slow proportional servo 
by themselves have already been tried and proven, the use of a servo- 
motor combining the two characteristics is suggested, the fast-acting 
flicker portion to alleviate quickly initial disturbances and to seciore 
a large response to an input signal. The main function of the propor- 
tional part is to secure stabilization around the neutral point of the 
range. 


The analysis contained herein pertains to one specific supersonic- 
model configuration for which satisfactory stability was achieved. It 
is very probable, however, that the optimum adjustments of the variable 
factors in the autopilot mechanism have not been realized, and it is 
believed that a more complete analysis utilizing an analog computing 
machine would show the optimum autopilot adjustment more closely. The 
results of the analysis contained herein show the effects of the fol- 
lowing conditions on the stability of the autopilot -model combination 
based mainly on the flight condition anticipated as a result of previous 
flight tests of the model: 




1 . Normal acceleration 
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2. Aerodynamic out-of-trim moment 

3. Static margin^ altitude^ and Mach number variation 

The analysis was continued further to include the response to an 
initial disturbance and the response to a command signal using approxi- 
mate physical autopilot transfer functions in the band of proportional 
autopilot operation. 


SYMBOLS 


6 canard control -surface deflection^ trailing edge down 

denotes positive deflection^ degrees 

6 error angle between gyroscope frame and instantaneous 

airframe attitude^ degrees 

K proportional-autopilot gain factor 

Q pitch angle measured from horizontal^ positive when nose 

is above the horizontal reference, degrees 

G first derivative of G with respect to time, degrees per 

second (d 0 /dt) 

a angle of attack, positive when nose is up relative to 

flight path, degrees 

a first derivative of a with respect to time, degrees per 

second (da/dt) 

7 flight -path angle, degrees (7 = 0 - a) 

g acceleration due to gravity, 32.2 feet per second per 

second 

n normal acceleration in g units 

M Mach mimber 

Y stability axis which passes through center of gravity and 

is perpendicular to vertical plane of symmetry 

ly moment of inertia about Y-axis, slug -feet square 
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L 

lift, pounds 

m 

mass, slugsj or pitching moment, foot pounds 

q 

dynamic pressure, pounds per square foot 

S 

wing area, square feet 

c,MAC 

mean aerodynamic chord, feet 

SM 

static margin 

V 

velocity, feet per second 

Cl 

lift coefficient (L/qS) 

Cjh 

pitching -moment coefficient (m/qSc) 

CL6 

SCl/S5 


SCl/^ 

Cm6 

^m/^S 


SCm/^ 

Cm^ 

acm/ae 


SCm/aa 

St 

value of control -surface deflection which co\interbalances 
out -of -trim moment, degrees 

mt 

out-of-trim moment caused by model misalinement , foot- 
pounds 

Cmt 

out -of -trim-moment coefficient (mt/qSc) 

t 

time, seconds 

T 

time-lag factor, seconds 

^ O 4- o ^ o 
^0 ^^5 ^^10 

time that the model attitude crosses the point in range 
denoted by subscript, seconds (for example, to° signi 
fies the instant of time that the model attitude 9 
is 0°) 
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o> 


c.p. 


C.g. 


D 

s 



frequency, radians per second 
model center of pressure 
model center of gravity 
differential operator (d/dt) 

Laplace -transform variable corresponding to differential 
operator 

signifies that inverse Laplace transform of function f(s) 
is to be taken 


Subscripts: 


input or forcing f^mction corresponding to a commajid 
calling for a change in attitude or to a sinusoidal 
input variation 

output or response fiinction, for example, system response 
to a command signal or to a sinusoidal input variation 


AIRFRAME AND AUTOPILOT PRINCIPLES AND DESCRIPTION 


A typical qualitative curve of the static servomotor displace- 
ment 6 plotted against the error signal e for a proportional -flicker 
autopilot system might be as follows: 
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For the theoretical analysis, the estimated pitch derivatives of 
the canard configuration, shown as a photograph an.d plan-view sketch in 
figure 1, were used. 

The autopilot system, as described in the following paragraph, 
constitutes one method of obtaining a proportional-flicker response, of 
which many variations are possible. It is also feasible that a com- 
pletely electric autopilot system could be devised which would function 
on the proportional-flicker principle. 

A schematic diagram of the type of automatic pilot investigated in 
this analysis is shown in figure 2. This system consists of a displace- 
ment gyroscope and a rate gyroscope which transmit error signals to a 
diaphragm by means of pneumatic Askanla pickoffs. (The signal is trans- 
mitted by an air jet to either of two holes in the pickoff block, which 
is connected to the servomotor diaphragm by rubber tubes.) The dia- 
phragm actuates the servomotor slide valve and the servomotor response 
becomes proportional to the input through the use of a feedback spring. 
The system thus described constitutes the proportional, part of the auto- 
pilot, which secures stabilization about the zero gyroscope reference 
point. The flicker portion of the system is obtained by use of elec- 
trical override pickoffs which are mounted on the displacement-gyroscope 
frame and which make contact at a preset angle of displacement with a 
pickoff attached to the outer gimbal of the displacement gyroscope. 

(The angle of displacement at which these pickoffs are set determines 
the switching point and the width of the band of proportional operation 
about the zero gyroscope reference point.) When either of the flicker 
pickoffs makes contact, one of the override solenoids is energized which, 
in turn, actuates the servo slide valve to cause the flicker action of 
the servomotor piston. The function of the leaf springs in the servo 
feedback linkage is to alleviate the feedback spring and diaphragm 
forces during flicker operation. With this arrangement, more rapid 
flicker action is assured with relatively small solenoids. In operation, 
a time lag occurs at the autopilot switching point. This lag causes a 
delay between proportional and flicker operation when the airframe atti- 
tude passes out of the proportional band and causes an overshoot under 
flicker operation when the attitude crosses the band from the other 
direction. It is possible to have delay and overshoot periods of differ- 
ent magnitudes. 


METHOD OF ANALYSIS 


The analysis contained herein consists of calculating the transient 
responses of the missile -autopilot combination for various initial condi- 
tions and approximating those conditions which might be encoiintered in 
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flight. The principle on which the proportional -flicker autopilot oper- 
ates is nonlinear; however, it does lend itself readily to aneilysis by 
means of linear approximations and a step-by-step solution. The first 
step is the response to a constant control deflection corresponding to 
flicker-servo action caused by an initial attitude disturbance outside of 
the proportional band. The airframe motion and corresponding autopilot 
response preceding the initial disturbance were not considered in this 
analysis. In the initial phase of the calculations, the autopilot oper- 
ation in the proportional band was approximated with a perfect con- 


trol 


5 =K. 

€ 


Calculations based on this assim^jtion are presented in the 


resiilts as approximations of the actioal autopilot -^odel pitch transient 
responses to an initial dist\jrbance . Since a physical servomotor will 
contain dynamic lag, further calculations were made using an approximate 


function j = f(D) for the autopilot response in the proportional band. 

Transient responses incorporating this approximation are also presented 
in the resiilts. 


Three methods of analysis were considered in obtaining the pitch 
transient responses of the supersonic model -autopilot combination pre- 
sented in the results. The solution for a constant control -surface 
deflection applies in the first step of each method. The first method 

can be generalized as the response to an Initial disturbance 
proportional band^ , the second as the response to an initial disturbance 
= f(D) in proportional bandj, and the third as the response to com- 
mand signal = f(D) in proportional band j . 


1 = K in 


Response to an Initial Disturbance 
^ = K in Proportional Bandj 


The forms of the equations of motion for constant speed and level 
flight are: 


57.3qSc 


D - CmQ D ^ 0 - + Cm^ d| a = Cmg 5 


57-3qS 


D 0 - ( „ D + CLcl la = Cls 5 


57.3qS 
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The general solution of these equations is broken down into two steps: 
( 1 ) 5 = Constant for the flicker portion of the solution and ( 2 ) & = Kg 
for the proportional part of the solution. For both steps the methods 
of Laplace (refs. 1 and 2) were used in obtaining the equation of 0 as 
a function of time^ 0 = f (t) • 


A qualitative example of a typical solution is as follows: 



band. Under this condition the control surface 6 is against its stop 
and the solution for 5 = constant applies until 6 crosses the pro- 
portional band. At this point a time-lag factor t is introduced to 
allow for the delays in the override solenoid and servomotor. After 
this overshoot period the control surface is assimied to move instantane- 
ously to a value which corresponds to 0 by the relation 6 = Kg = -K0. 
The control surface would actually move as shown by the dotted line. 
However^ the flicker response is estimated to be at a rate of at least 
700 ^ per second j therefore^ the step approximation is used because the 
solution is simplified^ as will be shown, without introducing an appre- 
ciable error. At the end of the delay period the proportional 
control-surface motion applies and the response is calculated according 
to the relation 6 = Kg. The initial conditions of this second step are 
obtained from the end conditions of the first step. Since the nature of 
the solution requires that the initial value of a be known, the tran- 
sient a = f(t) for a constant 6 was also derived in order to 
determine the value of a at the end of the first step. 


MCA TN 5^27 


9 


Response to an Initial Disturbance 
= f(D) in Proportional Band^ 

This method of solution was derived in order to define more com- 
pletely the motion of the 6 transient at the beginning of the second 
step. The solution for the first step is identical with that of the 
first method of analysis and the derivation for the solution of the 
second step is as follows. 

The relations between the airframe and autopilot parameters gov- 
erning the solution of the second step can be drawn in block diagram 
form as 


0^ 



Autopilot 

transfer 

function 

s 

Airframe 

transfer 

function 


0 

— ^ 






< 





where the transfer functions are functions of the differential operator. 

Using the methods of Laplace (refs, 1 and 2), the equations of motion 
are transformed to the form 


fl(s)0o(s) + f2(s)a(s) = Cmg5(8) + f3(s) 

^ 4 ( 3 ) 00 ( 2 ) + f 5 (s)a(s) = CLg5(s) + fg(s) 
where ^ 3 ( 3 ) and f 6 (s) contain the initial -condition terms. 

Expanding equations (1) and (2) and solving for 0o(s) yields 

fY(s)5(s) + fsCs) 


( 1 ) 

( 2 ) 


0n(s) = 


fo(s) 


( 3 ) 


It remains to define 5(s). By transforming the relation and using 

0i = 0 


6 = f(D)e = -f(D)0o 
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an equation of the form 

6(s) = 

is obtained where f]_i(s) contains the initial -condition terms 
involving 6 and . 

Combining equations (3) 9*^^ (^) solving for ^o(^) yields 


^o(s) = 




(5) 


where f]_2(s) contains the required initial values of 9, a, and 5 

and their derivatives. The inverse transform of equation (5) is the 
transient solution for the second step 


0(t) = l"^ 



By combining equations (3) and (U) and solving for 5(s), 
equation defining 6 as a function of s can be obtained^ the inverse 
transform of which is the second- step solution for the 6 transient 
response . 


Response to a Command Signal 
^ = f(D) in Proportional Band^ 

In an automatically stabilized missile, a co mm and signal can be 
obtained by changing the gyroscope reference point. This command signal 
may be generated, for example, by the radar unit of a target seeker, by 
an outside source in a guided missile, or by an altitude control. When 
the command signal is larger than the proportional band, the flicker 
portion of the autopilot under consideration will function, yielding a 
constant control-surface deflection, and the motion will be initiated 
toward the desired new attitude or new neutral point in the range. An 
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example of the Oq variation with time for this third method of 
analysis might he as follows: 



As is shown, the transient response is assumed to start from level 

flight and the new neutral point is outside of the original proportionaJ. 
hand. Initially the equations of motion are treated in the same manner 
as in the first step of the first and second method except that the 
initial value of Gq is zero. The transient response Qq to a 
constant 6 applies until Gq crosses the new proportional hand. At 
this point the time-lag factor t is again introduced to complete 
the first step. 

The method of the derivation for the solution of the second step 
found herein is carried through in general form in references 1 , 3 ^ 
and k. As mentioned previously, the autopilot frequency response is 
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approximated "by a transfer function of the form = f(D), where the 

error angle e is the angle between the gyroscope reference and the 
instantaneous airframe attitude. The second step^ shown as the dotted 
portion of the Go transient response, is the response to a step 9± 
where is the angle between the new neutral point and the value of Gq 

at the end of the first step, a new zero reference being taken for Gq 
at this point. The equations of motion are handled in the same manner 
as for the second step of the second method of analysis except that 
the relation governing the definition of 5(s) is 

6 = f(D)€ = f(D)(0i - 0o) 
therefore, equation ( 4 ) becomes 


5(s) = fiQ(s)|e^(s) - 

and the solution for 0o(s) takes the form 

fl4(s)ei(s) + fi2(s) 


0q(s) = 




( 6 ) 


The inverse transform of equation (6) gives the solution of the second 
step of the transient response to a command signal. 

The three aforementioned methods of analysis describe the general 
theories which were derived to obtain the solutions for the transient 
responses presented in this analysis. Any deviation from these general 
theories will be described in further detail in the results. 


RESULTS AND DISCUSSION 

Response to an Initial Disturbance 
= K for Second Step! 


Determination of the time-lag factor t .- The time-lag factor for 
the autopilot previously described is defined as the time that it takes 
for flicker-servo action after the flicker-actuating solenoid has been 
energized or deenergized. In order to obtain an approximate value of 
this factor for use in the calculations, a servomotor was constructed 
which operated on the proportional -flicker principle and which would be 
suitable for actuating the canard control surfaces of a model of the 
type shown as a photograph and sketch in figure 1. The results of an 
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experimental test of this servomotor are presented in figure 3- As is 
shown^ an average r of 0,03 second was obtained. It is not believed 
that more rapid flicker actuation can be obtained with this type of 
system due to the limitations on the size of solenoids j however^ it 
would not be difficult to increase the time-lag factor if necessary for 
stability. 

Determination of the proportional gain and amplitude of flicker- 
control deflection .- For these calculations the coefficients and model 
longitudinal pitch derivatives were estimated on the basis of sea-level 
flight at a Mach number of 1.8 with the static margin equal to 0.86c or 
approximately 14 inches. The geometric characteristics and estimated 
derivatives used are given in table I. 

A detailed description of the steps required for a solution such 
as shown in figure k is as follows: The first step is the response to 

a constant 5. The second step which starts at some time after the 
model attitude passes the proportional boundary is governed by the 
approximate relation 6 = -K0 for the major part of this analysis. 

This simplifies the calculations and still gives a good approximation 
of the actual 0 response. The values of the initial-condition terms 
necessary for the second-step solution of the equations of motion are 
obtained from their values at the end of the first step. If the 0 
response stays inside the proportional band and the motion dies out^ 
the solution is complete in two steps. However, if the model attitude 
again reaches the proportional band, the flicker-actuating solenoid will 
be energized, making more steps necessary. For the third step the 6 
response is assumed to hold a constant value which is determined from 
the relation 6 = -KO where Q is the proportional-band limit. The 
initial-condition terms are again found from their values at the end of 
the previous step. The time-lag factor determines the duration of the 
third step. When this period of time has elapsed the fourth step is 
initiated for which the constant flicker 5 applies as in the first 
step. This step-by-step solution is continued until the response of 
the model -autopilot combination is determined. 

The first calculated transient responses, presented in figure 
were based on a flicker b of ± 10 ^, a proportional band of ± 5 ^ which 
is designated in this and subsequent figures by the long-dashed lines, 
and a time-lag factor of 0.03 second. The initial pitch disturbance for 
the curves presented in this figure was 9^ and the value of K was -1 
for figures 4(a) and 4(b), while in figures 4(c) and 4(d) K = - 0 . 3 , 
where K is equal to the control- surface gain ratio b/e. Figure 4(a) 
shows the response when the first step ends at 0.03 second after the 
model attitude G passes the proportional boundary (+5^)^ designated 
by t^ + 0.03 second on the figure. Figure 4(b) shows the response 
when the first step ends at 0.03 second after the model attitude passes 
through the neutral point (0°) or tQ^ +0.03 second. This would 
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necessitate a more intricate flicker -actuating pickoff and relay system 
than that shown in figure 2> however, its development would not propose 
a major problem. As can be seen in both figures 4(a) and 4(b), the 
response is unstable as is indicated by the divergent transient 
oscillations. Curves for the same conditions are shown in figures 4(c) 
and 4(d) except that K = -0.3 "the proportional band. As can be 
seen, decreasing the gain in the proportional band does not increase 
the stability but actually has the effect of making the model-autopilot 
combination more unstable . 

Since it is apparent that a flicker control -surface gain of ±10° is 
too high for stability at M = 1.8, calculations for a lower flicker 
gain, namely ±5°, were made. Curves based on a flicker S of ±5° and 
on the same derivatives used for figure 4 are presented in figure 5. 

The curves of figure 5(a) show that, for an initial disturbance of 
Q = 9°, proportional band of ±5*^^ and for K = -0.3^ the transient 
response is undesirable when the first step ends at t^ + O.03 second, 
but when the switching point is extended to tQ° +0.03 second the tran- 
sient response is stable. The secondary oscillation which appears 
during the second step is induced by the instantaneous movement of the 
control surface 6. Figure 5(b) shows that, for a larger K (namely, 
for K = -l), the proportional band can be smaller (±3°) and the response 
will still be stable and damp out in two steps. For a smaller initial 
disturbance, however, the stability is more critical, as is shown by 
the curve which has an initial disturbance of ±5-5°* this reason 

the curve shown in figure 5(c) was calculated for an initial disturbance 
which is just outside of the proportional band, namely, an initial dis- 
turbance of 3.5° for a proportional band of ±5°- As is shown, the cal- 
culated transient response is divergent^ however, by increasing the 
time-lag factor t to 0.1 second as in figure 5(h) the transient is 
made to stay inside the proportional band in the second step. Actually, 
the autopilot in this instance could have two values of t, O .03 second 
when the attitude passes out of the proportional band and 0.1 second 
when the attitude crosses 0°. 

Since it appears that a stable transient response can be obtained 
by the proper adjustment of the vaxiable factors in the autopilot mech- 
anism, it was decided to increase the flicker control- surface gain 
to ±7° in order to obtain the advantage of a more rapid response time 
than that obtained with a gain of ±5°. For this reason the analysis of 
the calculated transient responses which follows will be based on a 
flicker 5 of ±7° and a proportional band of ±4°. The first responses 
for these conditions, shown in figure 6, were calculated using K = -1 
in the second step with the first step ending at some time after the 
model attitude has passed through the zero gyroscope reference. In this 
instance, the autopilot will function as a proportional system as long 
as the initial disturbance is within the +4° boundary. In the event that 
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the model attitude passes either side of these limits, the flicker 
control surface of ±7° will apply, whichever is corrective. The flicker- 
actuating solenoid will be deenergized and proportional control will 
again apply some time after the model attitude passes through the zero 
gyroscope reference. Proportional-flicker transient responses to a 9 
initial 6 disturbance were calculated with the first step ending at 
to° + 0.03 second and at tQ° + 0.1 second, as shown in figure 6. In 
this particular instance, the larger time lag yields a transient 
response which is closer to the zero 9 reference. The upper curve on 
figure 6 shows the response of a zero -phase -lag proportional autopilot 


ference between the response obtained from a linear autopilot and the 
responses obtained with the proportional -flicker autopilot. For the 
particular value of autopilot gain used in calculating the response 
obtained with the zero -phase -lag proportional autopilot, the 
proportional -flicker responses, although more oscillatory, are more 
rapid than the pure proportional autopilot response. The value of K 
used for the zero -phase -lag proportional autopilot response was - 0 . 6 . 

The choice of this autopilot constant was based on the method suggested 
in reference 4. This method is illustrated in figure 7 where the 
Nyquist diagram with a K of -0.69 is shown to be tangent to, and not 
greater than, the locus of points required to make the peak amplitude 
ratio of the model -autopilot combination have a value of 1.5* 

Acceleration effects .- Normal acceleration and 9, a, and 7 
transient responses for the model-proportional -flicker -autopilot com- 
bination are presented in figure 8 . The 9 transient is the curve with 
the first step ending at tQ° + 0.1 second, as shown previously in fig- 
ure 6 . As is shown, the trim a is about -2.5° and the steady-state 
normal acceleration is approximately -llg for the first step, with a 
peak acceleration overshoot of about -l 8 . 8 g. 

An indication of the model aerodynamic loads to be expected under 
flicker-autopilot operation can be obtained from the steady-state varia- 
tion of n /6 with static margin for a constant 6 input, as shown in 
figure 9 for sea -level flight. This curve is based on the variation of 
the model longitudinal derivatives with center-of -gravity location 
at M = 1.8. Figure 9 shows only the steady-state acceleration; however, 
the peak acceleration is the design factor, therefore acceleration 
curves such as the one shown in figure 8 are necessary to estimate the 
amount of acceleration overshoot. The method employed in deriving the 
relation on which the plot of n /6 against static margin is based is 
given in appendix A. 

The effect of out-of-trim moment .- The value of the out -of -trim 
moment is assumed to be that which will give a certain value of 6 
out of trim 8 ^ or 



This curve serves to illustrate the dif- 


l6 
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mt = CmgStqSc 


The out -of -trim-moment coefficient 



is therefore defined as 


Cmt - CiagSt 


This coefficient has the effect of adding another term to the moment 
equation as follows: 

(57^ ^ 

The solution of this equation, combined with the lift equation, 
for each step of the proportional-flicker transient responses shows the 
effect of introducing an assumed value of an out-of-trim moment. 

Transient responses, including the effect of an out -of -trim moment, 
are shown in figures 10 and 11. These figures are based on sea-level 
flight at M = 1.8 and with SM = 0.86c, and the initial 9 disturb- 
ance was 4.5 with K equal to -1 for the second step. It is shown in 
figure 10(a) that, without the out-of-trim effect, the mo del -autopilot 
transient response is made to damp out in the second step when the time- 
lag factor is increased to 0.1 second, that is, when the first step ends 
at 0.1 second after the model attitude has passed through the neutral 
point (0°). Therefore, the remaining curves shown in figures 10 and 11 
are based on the value to° +0.1 second for the length of the first 
step. Figure 10(b) shows the responses when out -of -trim -moment coeffi- 
cients of -0.012 and +0.012 are used. These values of Cmt cause the 
0 transients to trim about -0.57° and +0.57°^ respectively. However, 
the stability is not affected to any great extent since the response is 
damped to within approximately 1° of the trim value in 0.6 second in 
each case. More severe values of Cmt> namely, +O.O58 and -O.O58, were 
used in figure 11. These give trim 9 values of +2.87° and -2. 87°, 
respectively, which are only 1.13° from the proportional boundary 
limits (±4°). However, although the responses take longer to damp to 
the new trim values, they are still stable and stay inside the propor- 
tional band in the second step. For comparison, calculated transient 
responses including an out-of-trim moment and based on a zero-phase-lag 
proportional autopilot with K = -0.6 are also shown in figure 11. For 
these curves a Cmt °f ±0.035 was used because, for comparison, it is 
desirable that the trim 9 (±2.87°) "the same as for the proportional- 

flicker responses. In figure 11(a) it is shown that, since the 
proportional -flicker response must first pass through 0° and then 
oscillate back to the trim value, the response is not much faster than 
the pure proportional- autopilot response. However, in figirre 11(b) 
where the out-of-trim moment is negative, the proportional -flicker auto- 
pilot has the more rapid response time. 
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The response to an Initial disturbance Including Initial values 
of 9 and a . - An initial pitching velocity ("5^ of 50 degrees per 
second and an angle of attack (a) of 0.5° were used for the results 
presented in figure 12. These are representative of the values which 
might he encountered along the flight range of the type of model being 
used in this analysis. In figure 12(a) the initial 6 was 50 degrees 
per second with zero initial angle of attack^ while in figure 12(b) 
initial values of both 6 and a were used and the variation of normal 
acceleration and a with time are also shown. The stability of longi- 
tudinal transient responses is not affected greatly by including the 
effect of these initial conditions. 

The effect of Mach number change . - In order to determine the sta- 
bility characteristics of the model-proportlonal-flicker-autopilot com- 
bination for a different Mach number^ the derivatives for M = 1.4 
were estimated as given in table I for the same center -of -gravity loca- 
tion resulting in a static margin of 0.9c. Results based on these 
derivatives are presented in figure 13 as 6 and a transient responses 
to 9^ 4.5^ initial pitch disturbances. The 0 transients are com- 

parable^ except for Mach number^ to the curves with the first step 
ending at tQ^ + 0.1 second shown previously in figures 6 and 10(a) and^ 
although the response time is a little slower, the results of figure 13 
show stability. 

The effect of static-margin reduction . - The constants and aero- 
dynamic derivatives used for the curve presented in figure l4 were based 
on sea-level flight at M = 1.8 and SM = 0.3c and are given in 
table I. This amounts to a decrease in static margin from approximately 
l4 inches to 5 inches. As is shown, an undesirable result is obtained 
because the 0 transient response does not die out in the second step. 
Instead, the response diverges until what appears to be a pure flicker 
response is obtained. Therefore, in a system of this type, the problem 
of obtaining a more rapid response time is not solved by simply 
decreasing the static margin. 

The effect of altitude variation .- Except for altitude, the 
proportional -flicker response with the first step ending at tQ° +0.1 
second, shown in figure 6, is comparable to the curves presented in 
figure 15 where the pitch transient responses to a 9^ initial 9 dis- 
turbance are based on flight at 10,000 and 25,000 feet. As is shown, 
flight at altitude produces a slower response time but does not appre- 
ciably affect the stability of the model -autopilot combination because, 
in both cases presented in figure 15 , it is indicated that the 
transient response will die out in the second step. 
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Response to an Initial Disturbance 
= f(D) in Proportional Band^ 

Determination of the approximate autopilot functions . - The method 
employed to determine mathematical transfer functions which approximate 
the autopilot amplitude and phase response was as follows: First, the 

experimental autopilot amplitude and phase responses were obtained from 
oscillating-table tests of actual autopilots. Two autopilots were used 
in this analysis. The first consisted of the servomotor used previously 
to obtain the experimental results presented in figure 3^ and the 
oscillating-table error angle was measured by a German V-1 displacement 
and rate gyroscope which generates the input to the servomotor by the 
use of pneumatic Askania pickoffs. The autopilot amplitude and phase 
test points, thus obtained, were then plotted on semilog paper using a 
decibel scale for the amplitude response. Transparent templates based 
on plots of known quadratic and linear functions on this same semilog 
paper were then used until the combination of templates which most 
closely matched the autopilot test points was determined. The method 
used to make the log magnitude and phase templates and explanations of 
their applications are found in chapter 8 of reference 1. 

The approximate autopilot transfer function obtained by the afore- 
mentioned method is 


6 ^ 223(D + 27.2) 

^ + 141D + YJhh 

This function will then govern the proportional part of the proportional- 
flicker transient responses. A comparison of this approximate function 
with the autopilot test points based on the amplitude and phase response 
of the proportional-flicker servomotor obtained from oscillating-table 
tests of a V-1 displacement and rate gyroscope is presented in figure I6. 
The agreement between the approximate function and the experimental 
results can be seen in figure l6(a), where the approximate autopilot 
transfer function and the autopilot test points for oscillating-table 
amplitudes of ±1.13°, ±3-03°, and +4.85° are plotted on linear coordi- 
nates. The experimental results vary with oscillation amplitude due to 
the nonlinearities of the autopilot mechanism. However, the agreement 
between the mathematical function and the experimental results is sat- 
isfactory except for the low -amplitude oscillations (±1.13°) ^ where the 
amplitude response peaks more sharply. 

A further comparison between the approximate autopilot transfer 
function and the actual autopilot test points is made in figure l6(b) 
in the form of Nyquist diagrams, where the mathematical function and 
the test points are combined with the model transfer fimction to obtain 
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the loci of the polar plots. The model transfer fimction is based on 
sea-level flight at M = 1.8 and with SM = 0.86c. Satisfactory agree- 
ment is again shown between the approximate autopilot transfer function 
and the experimental results except for the Nyquist plot based on the 
±1.13° test points. 

The second approximate autopilot mathematical transfer function was 
based on the amplitude and phase response of a servomotor using a dis- 
placement gyroscope only to generate the input signal. This approximate 
function is 


5 ^ 2880 

€ + I56D + 3600 

and it is compared with the autopilot consisting of a displacement 
gyroscope plus servomotor in figure 17. Autopilot test points are 
plotted for oscillating-table amplitudes of ±4.2° and ±2.15 ) satis- 
factory agreement with the plot of the approximate mathematical function 
is shown. 

Autopilot containing displacement -plus-rate gyroscopes . - As shown 
previously in the method-of-analysis section, the function 5(s) and 
the initial values of the 5 derivatives are included in the derivation 
of the 0o(s) function in the second method of analysis. Although the 
manual solution for the 9 and S transient responses is much more 
involved when using this type of derivation, it is of value because it 
yields a more complete definition of the control -surface motion. 

Figure I8 shows the pitch transient responses to different initial 
9 disturbances based on the longitudinal derivatives given in table I 
for sea-level flight at M = 1.8 and with SM = 0.86c. The flicker 
gain is ±7° and the proportional band is ±4°. In figure l8(a), the 
initial 9 disturbance is 10° and the first step ends at 0.03 second 
after the model attitude has passed through the neutral point (0°). In 
this figure the calculated 6 response is also shown. In figure l8(b) 
the initial 9 disturbance is 4.5° with curves shown for the first step 
ending at tQ° +0.03 second and to° +0.1 second. The 6 response is 
plotted for the case where the first step ends at tQ° + 0.03 second. 

As is shown, the smaller time-lag factor (0.03 second) is more critical, 
with the second oscillation of the 9 transient actually crossing the 
proportional band. This would ordinarily necessitate handling the 
solution with more than two steps. However, it would be necessary to 
know the initial value of a for the third step and, since the charac- 

teristic equation for the a, transient when — = f(D) is tenth order, 

a Tnaminl solution was not attempted. It is more desirable to obtain a 
less oscillatory transient response which damps out in two steps by 
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increasing the time-lag factor, as is shown for the curve where 
T = 0.1 second. Figure l8(c) shows the 9 response to a 15^ initial 
disturbance, where for a t of 0.03 second the motion is damped in two 
steps . 

The agreement between the results shown for the second method of 
analysis using the autopilot containing displacement -plus -rate gyro- 
scopes and the results shown for the first method of analysis is good. 
The secondary oscillation which appeared in the second step of the first 
method of analysis, for example, is also present when using this second 
method of analysis, resulting in transient responses, the general 
shapes of which are compaxable. On this basis the validity of using 

the much simpler relation ^ = K for the major part of this analysis 
seems to be justified. 

Autopilot containing displacement gyroscope only .- The second step 
of the pitch transient responses to 10^ and 4.5^ initial 0 disturb- 
ances presented in figure I 9 is based on the transfer function 

6 ^ 2880 

^ + I56D + 3600 

These results show that stability can be obtained with a proportional- 
flicker autopilot made up of a servomotor and displacement gyroscope 
only. However, a comparison of figures I 8 and I 9 shows that the auto- 
pilot without rate yields a more oscillatory response in the 
proportional band. 


Response to a Command Signal 
= f(D) in Proportional Band^ 

The response to a command signal .- Pitch transient responses to 
4 . 5 ^ and, IQQ command slgna.ls are shown in figure 20. These responses 
are based on sea-level flight at M = 1.8 and with a static margin 
of 0.86c. The initial disturbance in each case is larger than the 
proportional band; therefore, the flicker-control setting of + 7 ^ applies 
for the first step of the solution, which ends at 0.1 second after the 
transient crosses the desired new attitude. In the second step the 
approximate transfer function 


6 

€ 


22^(D 4- 27 . 2 ) 

+ lUlD + 7744 
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applies. In figure 20(a) the desired change in attitude is 4.5^ or^ as 
explained in the method-of -analysis section^ the transient response 
shown is the response to a 4.5^ command signal. In figure 20(b)^ the 
command signal calls for a 10^ change in attitude. These results show 
that stable responses to command signals can be obtained with the 
proportional -flicker autopilot and the responses appear to be similar, 
except for direction, to the results for 4.5^ and 10 ^ initial disturb- 
ances shown in figure I 8 . A comparison between the proportional -flicker 
and the zero -phase -lag proportional response is made for both curves 
shown in figure 20 and in each case the proportional-flicker autopilot 
has the more rapid response time. 


CONCLUSIONS 


The proportional-flicker automatic pilot operates on a nonlinear 
principle, whereby a high-speed flicker servomotor response is combined 
with a low- speed proportional servomotor response for the purpose of 
obtaining stability and control in supersonic flight. Physcially, the 
autopilot motion operates about a zero reference within two bands. In 
the inner band, the autopilot output is proportional to the input and a 
flicker response overrides the proportional response at a fixed angle of 
gimbal displacement on either side of the zero gyroscope reference. The 
conclusions arrived at as a result of the analysis conducted herein, 
based on a specific supersonic model configuration, are as follows: 

1. Satisfactory stability can be obtained by the proper adjustment 
of the variable factors in the autopilot mechanism, namely, the propor- 
tional servo gain, the amplitude of flicker control deflection, the 
autopilot time-lag factor, and the point in the range where the autopilot 
switches from a flicker to a proportional system. 

2. A reasonable aerodynamic out-of-trim moment of the model will not 
affect the stabilization qualities of the proportional-flicker autopilot 
to any great extent. Decreasing the static margin appears to have more 

of an effect on the stability of the proportional-flicker — autopilot -model 
combination. For the same autopilot characteristics, a decrease in static 
margin from 0.86 to O .3 mean aerodynamic chord yields a transient 
response which does not die out but which diverges until what appears 
to be a pure flicker response is obtained. 

3 . Good agreement is shown between the response to an initial dis- 
turbance using a mathematical transfer function to approximate the auto- 
pilot and the method of assuming that the autopilot has a pure propor- 
tional response in the proportional band with instantaneous movement of 
the control surface at the switching point. 
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4. The proportional-flicker control system can be fabricated and 
appears to be a practical method for obtaining pitch stabilization of 
a supersonic pilotless aircraft. Therefore, trials of this system, 
paxticularly in supersonic vehicles, appear warremted. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., September I 5 , 1950* 
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APPENDIX A 

DETERMINATION OF THE RELATION WHICH EXPRESSES n/6 
AS A FUNCTION OF STATIC MARGIN 


The relation 


Cmg = Cls^ 


(x - SM) . de (y + SM) 

c d5 c 


( 1 ) 


expresses Cmg as a function of static margin. 


From the lift equation, the steady-state value of the normal 
acceleration in g can he derived for a constant control -surface 
deflection as 


n 


1 

32.2m 


-CLo,qS 


Cmg 

^ + CLgqS 


& 


where the steady state a/6 has been assumed to he equal to 
incorporating the relation 


^m5 

Cnkx’ 


By 


^ = SM 

Ct. c 
•^a 


the following relation is obtained, where the variation of the steady- 
state n/6 is expressed as a function of SM and Cmg 


n 


6 ■ 


_3S_ 

32.2m 






( 2 ) 


For a particular value of SM, Cmg can be evaluated from equa- 
tion (1). Then, by using these two values in equation (2), n/6 is 

determined. 
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Symbols not Previously Defined 


rate of change of canard lift coefficient with canard 
deflection^ 

^Law 

lift -curve slope of main wing^ hCjJ^Ou^ 

X 

distance from center of pressure of canard control surfaces 
to center of pressure of model, inches 

y 

distance from the center of pressure of main wing to center 
of pressure of model. Inches 

de/d6 

rate of change of downwash angle at wing due to deflection of 
canard control surfaces 

SM 

static margin, negative when c.p. is behind the c.g. 
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TABLE I 


ESTIMATED LONGITUDINAL DERIVATIVES 



derivatives in degree measure; 
m = L.922 slugs; c = 1.395 ft; 


ly = 37.66 slug-ft^; 

s = 2.52 


Static 

margin 

Mach 

number 

Cni0 


Cm5 


CL6 

CLa 

0.86c 

1.8 

-0.000138 

-0.052 

0.204 

-0.0000116 

-0.000051 

0.0607 

• 9 c 

1.4 

-.000187 

-.0621 

.023 

-.0000205 

-.00105 

.0692 

. 3 c 

1.8 

-.0000784 

-.0181 

.0182 

-.0000051 

-.000051 

.0607 


VARIATION OF FLIGHT CONDITIONS 


Altitude 

(ft) 

Mach 

number 

(lb/ft 2 ) 

V 

(ft/sec) 

sea level 

1.8 

4270 

1963 

sea level 

1.4 

2902 

1562 

10,000 

1.8 

3294 

1937 

25,000 

1.8 

1790 

1835 



(a) Photograph of model configuration. 

Figure 1.- Supersonic missile research model configuration. 
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(b) Plan -view sketch of model configuration. 
Figure 1.- Concluded. 
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Figure 2 .- Schematic diagram of a proportional- flicker autopilot system. 
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Figure 3 .- Proportional- flicker servo response to a sinusoidal input signal 

obtained from experimental tests. 
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(a) End of first step at 

t^o + 0.03 second; 

K = -1 for second step. 


(Id) End of first step at 

t^o + 0.03 second; 

K = -1 for second step. 


Figure 4.- Pitch transient responses to a 9 ° initial G disturbance based 
on sea-level fli^t at M = 1.8 and SM = 0.86c. Proportional band = ±5°^ 
T = 0.03 second, flicker 6 = ±10°. 


H 


VOVM 




(c) End of first step at 
t^o 0*03 second; 

K = -0.3 for second step. 

Figure 4.- 




(d) End of first step at 
tQ° + 0.03 second; 

K = -0.3 for second step. 
Concluded. 
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(a) Initial 0 disturbance = 9°> t = 0.03 second, proportional band = 

K = -0.3 for second step. 

Figure 5»- Pitch transient responses to an initial Q disturbance based 
on sea-level fli^t at M = 1.8 and SM = 0.86c. Flicker 6 = +5°. 
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(b) T = 0.03 second, proportional band = of first step at 

t^o + 0.03 second, K = -1 for second step. 

Figure 5«- Continued. 
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(c) Initial 9 disturteince = 3*5°^ t = O.O 3 second, 
proportional band = ±3°, end of first step at 
^0° 0*03 second, K = -1 for second step. 


(d) Initial 9 disturbance = 3*5'^^ t = 0*1 second, 
proportional band = ±3°^ end of first step at 
^0° ■*■ *^*^ second, K = -1 for second step. 


Figure 5* 
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Figure 6.- Pitch transient responses to a 9° initial 6 disturbance "based on 
sea-level flight at M = 1.8 and SM = 0.86c. Proportional "band = ±4°, 
flicker 6 = ± 1 °, K = -1 for second step. The calculated response of 
a zero-phase-lag proportional autopilot is shown for con 5 )arison. 
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Figure 7»“ Nyqulst diagram for a zero-phase-lag proportional— autopilot, model 
combination with K = -O .69 based on sea-level flight of the model and 
the estimated derivatives for M = 1.8 with SM = 0.86c. 



Figure 8 .- Normal acceleration n and longitudinal 0, a, eind 7 transient 
responses to a 9° initial 9 disturbance based on sea-level fli^t at 
M = 1.8 and SM = 0.86c. Proportional band = ±4°, flicker 6 = ±7°^ 
end of first step at tQO +0.1 second, K = -1 for second step. 
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Figure 9*- The steady-state sea-level variation of n/8 with static margin 
at M = 1.8 for a step 5 input based on the variation of the deriva- 
tives with center- of-gravity location. 
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(a) No out of trim. (b) End of first step at t^o + O-l second. 


Figure 10.- The effect of an out-of-trim moment on the pitch transient 
responses to a 4.5° initial B distiirbance based on sea-level flight 
at M = 1.8 and SM = 0.86c. Proportional band = -4°, flicker 
5 = ±7°, K = -1 for second step. 
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(a) Comparison of proportional- flicker autopilot with a zero-phase-lag 
proportional autopilot 0- = for a positive out-of-trim moment. 



(b) Comparison of proportional -flicker autopilot with a zero -phase -lag 
proportional autopilot 0 = K j for a negative out-of-trim moment. 

Figure 11.- The effect of an out-of-trim moment on the pitch transient 

responses to a initial 0 disturbance based on sea-level flight H 

at M = 1.8 and SM = 0.86c. Proportional band = flicker 

b = first step ends at tQO +0.1 second^ K = -1 for second 

stem . 
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(a) Initial 9 disturbance = 50 deg/sec; (b) Initial 0 distxirbance = 50 deg/sec; 
initial a disturbance = 0°. initial a disturbance = 0.5°. 

figure 12.- Longitudinal transient responses to a i^-.5° initial 9 disturbance including initial 

values of 9 and a based on sea-level flight at M = 1.8 and SM = 0.86c. Proportional 

band = ±4°^ flicker & = t7°^ first step ends at tQO +0.1 second, K = -1 for second step. 
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(a) Initial 6 disturbance = 9°* 



(b) Initial 9 disturbance = 


Figure 13. - Longitudinal transient responses to initial 6 disturbances 
of 9° QJid. 4.5° based on sea-level flight at M = 1.4 and SM = 0.9c. 
Proportional band = ±4°, flicker 6 = ±7°^ end of first step at 
tQO +0.1 second, K = -1 for second step. 
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Figure l4.- Pitch transient response to a 10° initial 6 disturb- 
ance based on sea-level fli^t at M = 1.8 and SM = 0.3c. 
Proportional band = flicker 6 = ±7°, end of first step 

at t^o + 0.05 second, K = -1 for second step. 
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(b) 25,000 feet. 

Figure 15.- Pitch transient responses to 9° initial 0 disturbances 
based on fli^t at 10,000 and 25,000 feet at M = 1.8 and 
SM = 0.86c. Proportional band = flicker 6 = ±7°, end of 

first step at t^o +0.1 second, K = -1 for second step. 
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(a) Arnplitude and phase response of the approximate autopilot function 
and of the actual autopilot test points. 

Figure 16.- Comparison of the function vhich approximates the proportional 
part of the proportional- flicker autopilot with the autopilot test 
points based on the an5)litude and phase response of a proportional- 
flicker servo obtained from oscillating-table tests of a V -1 displace- 
ment and rate gyroscope. 
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(b) Nyquist diagrams for approximate autopilot function-model combination 
and for longitudinal pitch oscillation amplitudes of tl.l3°, ^ 3 - 03 °, 
and ± 4 . 85 ° based on sea-level flight at M = 1.8 and SM = 0.86c. 


Figure I 6 .- Concluded 


Phc/se of/?g/ejC/eg A/T?p/ftude rat/o 


48 


NACA TN 3427 











Autopilot tost pO/MtS 
amplitude f dep) 

0 ± 4 .£ 

□ ± 2 .t 5 

Pu/?cf/oj7 H//j/ch approYimates autopilot 




























5 — 





Og 






















1 

G 

□ 
























f 

1 

























Figure 1 ?.- Con5)arlson of the function which approximates the proportional 
part of the proportional- flicker autopilot with the autopilot test 
points based on the amplitude and phase response of the servo obtained 
from oscillating-table tests of a displacement gyroscope only. 
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(a) Initial Q disturbance = lO^j end of first step at tQO + 0.03 second. 

Figure l8.- Pitch transient responses to an initial Q disturbance based 

on sea-level flight at M = 1.8 ajid SM = 0.86c. Proportional band = ±4^^ 

flicker 5 = y — = f(D) for second step including displacement-plus- ^ 

rate signal. ^ 
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(b) Initial 6 disturbance = 4 . 5 °* 
Figure l 8 .- Continued. 
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Figure l8.- Concluded. 
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Figure 19.- Pitch transient responses to 10° and 
B disturbances based on sea-level fli^t at 
SM = 0.86c. Proportional band = flicker 

of first step at t^o +0.1 second, j = f(D) 
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M = 1.8 with 
6 = ±7°, end 


step using displacement signal only. 


for second 


NACA TN 5427 


NACA-Langley - 5-17-55 - 1000 


i 


J » 


f 






Figure 20.- Pitch transient responses to and 10° command signals 

based on sea-level flight at M = 1.8 with SM = 0.86c. 
Proportional band = ±l4-°, flicker 6 = ±T°» ^ = f“(P) for second 

step including displacement-plus- rate signal. The responses of a 
zero-phase- lag proportional autopilot are shown for con 5 )arison. 
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